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ABSTRACT 

The nuclei of active comets emit molecules anisotropically from discrete vents. As the 
nucleus rotates, we expect to observe periodic variability in the molecular emission line 
profiles, which can be studied through mm/submm spectroscopy. Using this technique 



Based on observations carried out with the IRAM 30-m, JCMT 15-rn, and CSO 10.4-m telescopes. IRAM is 
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supported through NSF grant AST-0540882 (USA). 
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we investigated the HCN atmosphere of comet 103P/Hartley 2, the target of NASA's 
EPOXI mission, which had an exceptionahy favorable apparition in late 2010. We 
detected short-term evolution of the spectral line profile, which was stimulated by the 
nucleus rotation, and which provides evidence for rapid deceleration and excitation of 
the rotation state. The measured rate of change in the rotation period is +1.00 it 0.15 
min per day and the period itself is 18.32 it 0.03 hr, both applicable at the epoch of 
the EPOXI encounter. Surprisingly, the spin-down efficiency is lower by two orders 
of magnitude than the measurement in comet 9P/Tempel 1 and the best theoretical 
prediction. This secures rotational stability of the comet's nucleus during the next few 
returns, although we anticipate a catastrophic disruption from spin-up as its ultimate 
fate. 

Subject headings: comets: general — comets: individual (103P) — radio lines: general 



Introduction 



103P/Hartley 2 (hereafter 103P) is a Jupiter Family comet (JFC), which currently has a 
6.47-year orbital period and perihelion at 1.06 AU. On UT 2010 Oct. 20.7 it reached the mini- 
mum geocentric distance of only 0.12 AU, making by far the closest approach to the Earth since its 
discovery ( Hartley|1986 ), and becoming a naked-eye object. Shortly after, on UT 2010 Nov. 4.6, the 
comet was visited by NASA's EPOXI spacecraft, which approached the nucleus within 700 km and 
provided detailed images and spectra ( A'Hearn et al.||2011 ). Both the ground-based data, obtained 
at the unusually favorable geometry, and the unique observations carried out by the spacecraft, 
create an exceptional platform for new groundbreaking investigations. 

Taking advantage of this unique opportunity we addressed one of the burning problems of 
planetary science, which is the characteristic lifetime of active cometary nuclei. They can be 
dynamically ejected out of the sublimation zone, or decay through collision, tidal disruption, de- 
volatilization, deactivation, disintegration, and - the least understood - spontaneous fragmentation 
( |Jewitt|[l992l ). Of several scenarios proposed to explain the last process, a notable one is rotational 
break-up, which must occur for spinning bodies when the centripetal force surpasses self-gravity 



and material strength (e.g. Davidsson 2001). Interestingly, models show that torques exerted by 



gas jets can significantly increase or decrease the rotation rate during a single orbit (e.g. Gutierrez 



et al. 2002), although the magnitude strongly depends on the nucleus size and activity, and the 



poorly understood effective moment arm. The first unambiguous measurement of this effect was 



obtained only recently for comet 9P/Tempel 1, which was found to be slowly spinning up (Belton & 
Drahus|2007 ; Belton et al.|20lT ) . Clearly, further observational studies of the dynamics of cometary 
rotation are essential to establish the nucleus lifetime with respect to rotational break-up. 

Since the influence of jets must be strongest for small, elongated, and active comets (e.g. 
Samarasinha et al. 1986, Jewitt|[l999 Drahus & Waniak 2006), the nucleus of 103P is an excellent 
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test body: it is one of the smallest known JFCs (equivalent radius ~ 600 m), it is elongated and 
active ( Lisse et al.||2009 A'Hearn et al.|[20lT ), and hence was predicted to show measurable period 
changes during the last return (Drahus 2009 Lisse et al. 2009). In general, however, it is very 
difficult to investigate the rotation of an active nucleus; especially of a very small one and from 
the ground - securing simultaneously the necessary measurement precision. Our earlier studies 



show that perhaps the best possibilities are offered by mm/submm spectroscopy (Drahus et al. 



2010). When the molecules are emitted from discrete vents, or the body is irregular, the diurnal 



rotation modulates the emission line profiles, which can be conveniently studied in a velocity- 
resolved spectral time series. The effect will be best visible when every single spectrum results from 
the molecules released at similar nucleus rotation phases, which is satisfied when the molecules leave 
the projected telescope beam on a timescale much shorter than the body rotation period. The best 
candidates are hence nearby comets rotating slowly, and observed ideally in small beams using short 
integration times with the largest antennas. Consequently, we identified 103P as a terrific target for 
large ground-based mm/submm telescopes ( Drahus||2009 ), owing to the unusually small geocentric 
distance, localized activity, and more than sufficiently slow rotation (the last two confirmed after 
the fact). 



2. Observations 

Between early September and mid December 2010 we used all large single-dish ground-based 
mm/submm facilities operating around one millimeter and offering open time, to carry out spectral 



monitoring of HCN. This molecule is a particularly good tracer of the cometary rotation (Drahus et 
al.|2010 ) thanks to its well-established origin directly from the nucleus and because it has by far the 
brightest emission lines in the one-millimeter atmospheric windows. Our campaign was the first 
at these wavelengths specifically designed to investigate the rotation state of a comet, and thus 
provided an unprecedentedly rich and dense velocity-resolved spectral time series of a cometary 
parent molecule. 

In this work we use 438 spectra of the J(3-2) and J(4-3) rotational transitions (Fig. [T]), 
collected with IRAM, JCMT, and CSO on 20 nights between UT 2010 Sep. 29.3 and Dec. 15.6 
(Table [T]). The spectra were calibrated in the standard manner; each covers 10 to 15 min, and 
was analyzed at spectral resolutions of 0.10 km s'^ (high resolution; HR) and 0.25 km s^^ (low 
resolution; LR). Details of the observations and reductions will be presented in a subsequent paper. 



3. Analysis 

The HR spectra were first parameterized in terms of the line area and center velocity and then 
converted to physical quantities: HCN production rate Q and median radial HCN gas-flow velocity 
Urad- All these steps were performed following our standard methods and procedures ( Drahus|2009 
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Drahus et al. 2010), and will be described in a subsequent paper. In this work we only highlight 



our approach to estimate errors (Appendix [A|) . 

Although unavoidably affected by model simplifications, the physical quantities are naturally 
free (to the first order) of trivial instrumental, topocentric, and transition-specific differences, hence 
making the data set fairly homogeneous and suitable for the investigation of variability. Figure [2] 
shows the evolution of Q with heliocentric distance r, and displays frad against phase angle (j). 



The first trend in Fig. [2] was measured for several comets before (e.g. Biver et al. 2002a 



Drahus et al.|2010D . We see an increase in the production rate as 103P approaches perihelion and a 
decrease afterwards, with superimposed short-term variability responsible for the excessive scatter. 
The maximum mean-diurnal level of HCN is about 1 x 10^^ molec s~^, which seems high for a 
JFC but not abnormal for a comet (Fig. [s]). The ratio of HCN to H2O ( [Biver et al]|2010t [A'Hearn 



et al. 2011) at the level of 0.1% is also typical of comets (Biver et al. 2002b, Bockelee-Morvan et 



al. 2004). After subtracting the heliocentric trend we derived production-rate deviations Alog(Q) 
about the mean-diurnal level. 

The expected correlation of Vi-ad with (f) has been observationally confirmed only recently 



(Drahus et al. 2010); however, in our data we do not find any trend (Fig. [2]), because of the 
highly limited phase-angle coverage, short-term variations, and observing noise. 



4. Nucleus Rotation 



While the short-term variations could be sporadic ( A'Hearn et al.|2005 ), it is also entirely pos- 



sible that they were stimulated by the nucleus rotation. We pursued this issue using an innovative 
technique: the Dynamized Structural Periodicity Analysis (DSPA), which we created specifically for 
this purpose (Appendix [B]) . The method was developed from our Dynamized Periodicity Analysis 
(DPA), which - utilizing classical algorithms for a time series of signals - introduces a dynam- 



ical formula for the body rotation and rigorous weighting of the input data (Drahus &; Waniak 



2006). DPA determines the rotation frequency simultaneously with its time-derivative from the 
entire data set, taking advantage of the maximum available time base, and therefore is far superior 
compared to separate classical (constant) period determinations from smaller subsets. Recently 
DPA demonstrated excellent sensitivity in providing important evidence for the slow spin-up of 
comet 9P/Tempel 1 ( Belton &: Drahus||2007 Belton et al.pOlT ). DSPA goes further, and analyzes 



periodicity not only in the level, but in the complete structure of the signal in the input data, prop- 
erly combining the information from independent data channels. At the same time it keeps all the 
properties of DPA, and for a single-channel time series fully reduces to the latter. The technique 
does not rely on any variability models, neither versus time nor across the channels. Although 
DSPA allows for any a-priori law controlling the effective torque, in this analysis we assumed it was 
constant. 



Using the derived production-rate variations Alog(Q), we normalized the original spectra to 
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eliminate all the trivial effects and remove the physical heliocentric trend (Fig. [2]). In order to 
secure an affordable computation time and reasonable S/N in the individual data channels, DSPA 
was applied to the normalized LR spectra. We used 17 channels from the velocity interval between 
-2 and +2 km s~^, which is more than enough to adequately sample the line profile. They were 
analyzed with two well-established kernels: weighted version of Phase Dispersion Minimization 
(PDM) and weighted Harmonic Fit (HF) ( [Drahus k Waniak|[2006 ) . 



Figure [4] presents dynamical periodograms resulting from the analysis. The best single-peak 
solution is found for the rotation frequency /o = 0.0546 it 0.0001 h~^ (given at the epoch of 
the EPOXI encounter; indicated by index zero) and frequency time-derivative d//dt = (—2.0 it 
0.3) X 10-6 h-2. It is equivalent to the period of 18.32 ± 0.03 h increasing at the relative rate of 
0.07% ±0.01% (e.g., 1.00 ±0.15 min/day, or 0.75 ±0.10 min/cycle), both applicable at the EPOXI 
encounter. The result is insensitive to the uncertainties resulting from the earlier reductions - 
including the removal of the heliocentric trend, as the data most affected enter the analysis with 
very low relative statistical weights. 

We note that the global solution is very close to the third multiple of the above single-peak 
solution, which means that the pattern of mass loss repeats significantly better every three rotation 
cycles. We explain this behavior by excitation of the rotation state: every three cycles of the 
fundamental periodicity the independent rotation modes coincide at similar phases and then the 
3-cycle repeats. But the resonance is not perfect. The spectra from the same rotation phases 
and corresponding cycles still show differences in the line profiles (e.g.. Fig. [T] shows spectra from 
phase 0.5 and Cycle B [cf. Fig. |4j, corresponding to the moment of the EPOXI encounter). We 
attribute this effect primarily to small differences in the instantaneous nucleus orientation, which 
cause differences in the direction of gas flow. Interestingly, the larger nucleus of comet 9P/Tempel 1 
was found close to a fully- relaxed rotation state ( Belton et al.||20lT ), while the even bigger nucleus 



of IP/Halley was in a significantly excited state (Belton et al. 1991). 



Figure [4] also shows the production-rate variations Alog((5) and median radial gas-flow ve- 
locities t>rad phased according to the best single-peak solution and its third multiple. The line 
variability comes primarily from the diurnal changes in the production rate, which follow a nearly 
sinusoidal trend. The single-peak solution is consistent with a single strong active area being the 
primary source of HCN. As the EPOXI encounter occurred at the rotation phase of 0.5, when the 
HCN production rate was shortly before the diurnal maximum, we suppose that this area coincides 



with the strongest jet observed by the spacecraft (A'Hearn et al. 2011). This also implies that the 
higher-order multi-peak solutions, including the triple-peak third multiple indicated by the global 
minimum, cannot be plausibly associated with the fundamental periodicity. 
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5. Discussion and Conclusions 



Our best solution indicates a rapid spin-down of lOSP's nucleus. Spin-down at a comparable 



rate was suggested earlier for comet C/2001 K5 (Drahus & Waniak 2006), but the only robust 



measurement, obtained for 9P/Tempel 1, indicates spin-up at a 20-times smaller rate (Belton & 



Drahus 2007; Belton et al. 2011). Dynamics of the rotation rate can be interpreted in terms of a 



simple model (e.g. Samarasinha et al. 1986; Jewitt 1999; Drahus Sz Waniak 2006), where the rate 



of change in frequency d//dt is given by: 



df/dt 



15 Vs Q 



tot 



16^2 ^4 ^ 



(1) 



where Vs is the gas sublimation velocity, Qtot is the (mean-diurnal) production rate of all gaseous 
species, Ry is the volume-equivalent nucleus radius, p is the nucleus bulk density, and k is the 
dimensionless effective moment arm which is a measure of acceleration efficiency; k is negative 



for spin-down and positive for spin- up, and may exceed unity for non-spherical bodies (Drahus & 

6 h-2, and Qtot = QHa O+QcOa = 450 
^ (A'Hearn et al. 



Waniak 



2006). Substituting the measured d//dt = -2.0 x 10" 



kg s"^ dA'Hearn et al.|[20Tl| , = 400 m s^^ ( |Combi et aLl|2004[ ), 



;tot 

400 kg m^ 



2005] ), and Ry = 580 m ( |Lisse et al. | |2009| [A'Hearn et"aL][20TT| ), we find k = -0.04%. The same 
equation yields k = +4% for the properties of comet 9P/Tempel 1 (i?v = 3000 m [A'Hearn et al. 
2005], Qtot = Qh20 + Qcoa = 160 kg s"^ [Feaga et al. 2007], and df/dt = +9.7 x 10"^ h'^ 
[Belton et al. 2011]; otherwise the same as above), which is consistent with a model prediction 
\k\ ~ 5% ( |Jewitt|p99l ). Even though K in 103P is two orders of magnitude lower compared to 9P, 



the comet changes its rotation frequency a factor of 20 faster. This is understandable from Eq. ([T]), 
since df/dt has a strong dependence on the nucleus size, and the nucleus of 103P is tiny compared 
to 9P. Moreover, 103P is more active per unit surface area (Fig. |3]) and approaches the Sun closer, 
therefore it has a higher total sublimation rate at perihelion, in spite of the much smaller total 
area. 

If the pattern of mass loss remains stable, the above result lets us crudely speculate about the 
rotational past and future of 103P on the current orbit. We estimate the orbit-integrated frequency 
change at the level of 0.012 h~^ and assume the critical rotation period of approximately 3 h 
( Davidsson|2001 ) . This allows for no more than ~ 20 previous revolutions. At the present rate, the 
nucleus will stop rotating during the fourth or fifth return from now (perihelion in 2036 or 2043). 
If it survives the increased thermal stress, it should start spinning up in the opposite direction, 
and reach the rotational disruption limit ~ 25 orbits later (year ~ 2200). We note that rotational 



break-up has by far the shortest timescale determining the lifetime of such a small object (Jewitt 



1992, 1999 Lisse et al. 2009). Moreover, bearing in mind the location of the main active area on 



the nucleus (A'Hearn et al. 2011), even a small change in the jet configuration may dramatically 
increase the effective moment arm - bringing 103P to the end of its life within one or few orbits. 



We expect that small, young, and naturally volatile-rich comets, having typical effective mo- 
ment arms of the order of 5%, routinely experience rotational instability and fragmentation right 
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after being injected to the sublimation zone yet before being discovered. In this way the observed 



flat size distribution of comet nuclei (e.g. Snodgrass et al. 2011) can be naturally explained. With 



its tiny and active nucleus, 103P appears as a lucky survivor, protected for some limited time from 
rotational disruption by the remarkably low factor k. 

We thank M. Polihska for providing resources and the telescope staffs for their excellent work. 
This project was supported by NASA through the Planetary Astronomy Program grant to D. Je- 
witt. 

Facilities: IRAM 30-m (EMIR), JCMT (HARP), CSO (Z-Rex). 



A. Error Evaluation 

We have examined all error sources and concluded that the signal uncertainty in the individual 
spectral channels is dominated by observing noise and imperfect pointing. 

Observing noise was estimated from signal variations about the baseline assuming a normal 
distribution. Since noise is nearly entirely controlled by the sky and receiver temperatures (sur- 
passing the comet antenna temperature by about two orders of magnitude), we assumed the same 
noise level inside the line window. Noise dependence on frequency is negligible in this extremely 
narrow spectral range. 

The influence of pointing errors was estimated from a simple model, which we have developed 



for this purpose from our earlier constructions (Drahus et al. 2010). For each telescope axis we 



assume a normally-distributed pointing error with a standard deviation of 2" for IRAM and JCMT, 
and 3" for CSO; it is easy to show that this is equal to the most probable total pointing offset. We 
also adopt a simple steady-state isotropic coma model with inflnite molecular lifetime and brightness 
proportional to the number of molecules. Under these assumptions we calculated Probability 
Density Functions (PDFs) of signal loss for each telescope, and neglected other (less signiflcant) 
effects from imperfect pointing. 

To calculate the errors of all the parameters of interest we followed our Monte-Carlo approach 



(Drahus & Waniak 2006 Drahus et al. 2010), and generated 500 simulations using the derived PDFs. 
A starting point for each simulation was an original spectrum corrected for the most probable signal 
loss from pointing. To each spectral channel we added a random realization of noise and then scaled 
the complete spectrum according to a random realization of signal loss from pointing. Then, by 
analyzing the parameters' variations about their original values, we derived their standard errors. 

The signal errors in the individual spectral channels feature several interesting properties. 
With signal approaching zero they asymptotically approach the noise limit, and oppositely, they 
approach the pointing limit with the signal increasing to inflnity. The pointing component is 
proportional to signal whereas the noise component does not depend on signal. While the noise 
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component can be controlled (reduced) by spectral or temporal binning, the pointing component 
cannot, because it affects all the channels in one spectrum in the same way, and similarly several 
spectra taken within short time (the latter not being taken into account by our simple model) . For 
this reason the pointing component should not be included when analyzing signal variability across 
the channels. Moreover, the error from pointing is nearly one-sided (positive) while the one from 
noise is symmetric. 

Because several of the assumptions of our model may not be well satisfied, the provided uncer- 
tainties should be considered as crude estimates. Nevertheless, to our best knowledge it is the first 
approach to quantify how the imperfect pointing of radio telescopes affects the observed spectra. 



B. Dynamized Structural Periodicity Analysis 

The DSPA method is applied to a time series of M measurements. Quality of data phasing 
is measured by a function 0, which is a variance ratio of the phased and unphased data. DSPA 
calculates functions 9i{f,df/dt) of DPA (Drahus & Waniak 2006| > for every data channel i, and 



averages the results from the considered N channels to provide 9{f,df/dt): 

N 



N 
i=l 



(Bl) 



using weights Wi defined as square-mean signal over mean-square error: 



a? 




W^ = ^={'-^^ 1 (B2) 



where Sij is the signal in i-th data channel and j-th time moment, aij is its standard error, and 
Wij = o^j' ■ Note that both functions Bi used by DPA, i.e. the weighted version of PDM and the 
weighted HF, use the same weights Wij as above. Since our signal errors are asymmetric, we used 
their mean-square values for this purpose. 

Weighting with Wi implemented in DSPA naturally prefers channels with strong signal, making 
the method insensitive, from the theoretical point of view, to the selection of line window (as long 
as the line is fully inside); channels having on average zero signal contribute to Q with zero weights. 
Moreover, it prefers channels with small uncertainties, which is particularly important in this 
application since imperfect pointing differentiates signal error across the channels; for identical 
uncertainties in every channel this term cancels out. It is therefore clear that in our case the 
effective channel weight is established by the fine balance between channel signal and noise - as 
the error from pointing is proportional to the signal. 



-9- 



REFERENCES 

A'Hearn, M. F., et al, 2005, Science, 310, 258 
A'Hearn, M. F., et al., 2011, Science, submitted 
Belton, M. J. & Drahus, M., 2007, BAAS, 39, 498 

Belton, M. J. S., Mueller, B. E. A., Julian, W. H., Anderson, A. J., 1991, Icarus, 93, 183 
Belton, M. J. S, et al., 2011, Icarus, in press 
Biver, N., et al., 2002a, EM&P, 90, 5 
Biver, N., et al., 2002b, EM&P, 90, 323 

Biver, N., Bockelcc-Morvan, D., Crovisier, J., Lecacheux, A., Frisk, U., Floren, H.-G., Hjalmarson, 
A., 2010, CBET, 2524 

Bockelee-Morvan, D., Crovisier, J., Mumma, M. J., Weaver, H. A., 2004, in Comets II, eds. M. C. 
Festou, H. U. Keller, H. A. Weaver (Tucson: Univ. of Arizona Press), 391 

Combi, M. R., Harris, W. M., Smyth, W. H., 2004, in Comets II, eds. M. C. Festou, H. U. Keller, 
H. A. Weaver (Tucson: Univ. of Arizona Press), 523 

Davidsson, B. J. R., 2001, Icarus, 149, 375 

Drahus, M., 2009, Microwave observations and modeling of the molecular coma in comets, Ph. D. 
thesis (Univ. of Gottingen) 

Drahus, M. & Waniak, W., 2006, Icarus, 185, 544 

Drahus, M., Kiippers, M., Jarchow, C, Paganini, L., Hartogh, P., Villanueva, G. L., 2010, A&A, 
510, A55 

Feaga, L. M., A'Hearn, M. F., Sunshine, J. M., Groussin, O., Farnham, T. L., 2007, Icarus, 191, 
134 

Gutierrez, P. J., Ortiz, J. L., Rodrigo, R., Lopez-Moreno, J. J., Jorda, L., 2002, EM&P, 90, 239 
Hartley, M., 1986, lAU Circ. 4197 

Jewitt, D. C, 1992, in Proceedings of the 30th Liege International Astrophysical Colloquium, eds. 
A. Brahic, J.-C. Gereard, J. Surdej (Liege: Univ. Liege Press), 85 

Jewitt, D., 1999, EM&P, 79, 35 

Lisse, C. M., 2009, PASP, 121, 968 



-10- 



Richardson, J. E., Melosh, H. J, Lisse, C. M., Carcich, B., 2007, Icarus, 190, 357 

Samarasinha, N. H., A'Hearn, M. F., Hoban, S., Klinglesmith, D. A., 1986, in ESA Proc. of the 
20th ESLAB Symp. on the Exploration of Halley's Comet (ESA SP-250, vol. 1), 487 

Snodgrass, C, Fitzsimmons, A., Lowry, S. C, Weissman, P., 2011, MNRAS, in press 



This preprint was prepared with the AAS IM^jX macros v5.2. 



Table 1. Journal of observations 



UT Date 2010 


Telescope Cover.'* Num.'' r'^ 


Ad 


(j)^ Beam^ Escape^ 


(mid time) 


(h) (AU) 


(AU) 


(°) (km) (mill) 



Sep. 


29.4795 


JCMT 


6.6 


20 


1 


.1310 


0.1926 


44.1 


958 


20.0 


Sep. 


30.4229 


JCMT 


8.1 


28 


1 


.1265 


0.1873 


44.5 


932 


19.4 


Oct. 


1.2865 


JCMT 


2.8 


11 


1, 


.1225 


0.1826 


44.8 


908 


18.9 


Oct. 


16.6228 


JCMT 


6.5 


21 


1, 


.0710 


0.1241 


50.6 


617 


12.9 


Oct. 


17.5706 


JCMT 


10.0 


33 


1, 


.0691 


0.1228 


51.1 


611 


12.7 


Oct. 


18.5236 


JCMT 


5.0 


16 


1, 


.0673 


0.1218 


51.6 


606 


12.6 


Oct. 


19.5997 


JCMT 


8.3 


27 


1 


.0655 


0.1211 


52.3 


602 


12.5 


Oct. 


23.6593 


CSO 


6.7 


20 


1, 


.0606 


0.1225 


54.7 


1249 


26.0 


Oct. 


24.5773 


CSO 


8.3 


24 


1, 


.0599 


0.1236 


55.2 


1260 


26.3 


Nov. 


1.6117 


JCMT 


6.1 


20 


1, 


.0604 


0.1447 


58.5 


720 


15.0 


Nov. 


2.2058 


IRAM 


7.9 


24 


1, 


.0609 


0.1469 


58.6 


470 


9.8 


Nov. 


3.1957 


IRAM 


8.5 


26 


1 


.0619 


0.1506 


58.7 


482 


10.1 


Nov. 


4.1974 


IRAM 


8.4 


25 


1 


.0631 


0.1546 


58.8 


495 


10.3 


Nov. 


5.1820 


IRAM 


8.7 


28 


1 


.0644 


0.1586 


58.8 


508 


10.6 


Nov. 


10.6112 


CSO 


7.8 


24 


1, 


.0748 


0.1829 


58.1 


1864 


38.8 


Nov. 


11.5986 


CSO 


8.2 


30 


1, 


.0773 


0.1876 


57.8 


1912 


39.8 


Nov. 


12.6163 


CSO 


7.8 


28 


1, 


.0800 


0.1924 


57.5 


1962 


40.9 


Nov. 


13.5404 


JCMT 


0.8 


4 


1, 


.0826 


0.1969 


57.1 


980 


20.4 


Nov. 


13.5903 


CSO 


7.4 


25 


1 


.0827 


0.1972 


57.1 


2010 


41.9 


Dec. 


15.6115 


JCMT 


0.8 


4 


1, 


.2471 


0.3650 


37.8 


1816 


37.8 



Note. — IRAM and CSO observed HCN J(3-2) and JCMT observed HCN J(4-3). 

''Time coverage between the first and last spectrum. 

''Number of collected spectra. 

"^Heliocentric distance. 

''Topocentric distance. 

^Topocentric phase angle. 

'HWHM of the main beam at comet distance; at the observed frequencies the HWHM is 4.4" for 
IRAM, 6.9" for JCMT, and 14.1" for CSO. 

^Minimum escape time from the main beam, needed to reach HWHM with the assumed constant 
velocity of 0.8 km s~'. 
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Fig. 1.— 



Example of HR spectra. The unshaded region was used for calculations. 
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Fig. 2. — Orbital trends in the data. Top panel: Dependence of the production rate Q on heho- 
centric distance r (normahzed to the perihehon distance q) in a log-log space. Negative log(r/g) 
indicate positions before perihelion. The solid line shows a semi-empirical fit that was used to re- 
move the heliocentric trend and calculate the production-rate deviations Alog((5). Bottom panel: 
Behavior of the median radial velocity Ui-ad with phase angle (p. Note: The EPOXI encounter 
occurred at log(r/g) = 0.002 and (f) = 58.8°. 
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Fig. 3. — Compilation of the normalized HCN production rates Qiau for comets with known area- 



equivalent nucleus radius Ra, presented in a log-log space (see Drahus et al. 2010, and references 
therein). The production rates were normalized to the common heliocentric distance r = 1 AU 
using the canonical Q{r) ~ whenever needed. The solid line is a square function fitted in a 
log-log space. It represents a typical (constant) active fraction of the nucleus. 
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Fig. 4. — Results of DSPA. Upper panels: Dynamical periodograms for the regions of investigated 
space. Linear color scale denotes ^disp 

= log{r]{6 — 6*111111 )/(^max - 6*111111) + 10 — ry), where r/ = 9 
and function 9 = 9(f,df/dt) is defined in Appendix [b] Lower panels: Rotation phase profiles for 
Alog(Q) and Urad for two periodicity solutions (see text). The shaded regions indicate duplications 
of the phase space. The EPOXI encounter occurred during Cycle B at the rotation phase of 0.5. 



